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We present a mathematical framework for simulation of optical fields in complex gravitational-wave 
interferometers. The simulation framework uses the two-photon formalism for optical fields and 
includes radiation pressure effects, an important addition required for simulating signal and noise 
fields in next-generation interferometers with high circulating power. We present a comparison of 
results from the simulation with analytical calculation and show that accurate agreement is achieved. 



I. INTRODUCTION 

Next-generation gravitational-wave (GW) interferom- 
eters, such as those planned for Advanced LIGO Q, are 
designed to have a fifteen-fold improvement in sensitiv- 
ity over present-day detectors Q . Among the techniques 
planned to achieve this improved sensitivity is an increase 
in the input laser power. The higher laser power re- 
duces the shot noise limit at frequencies above ~100 Hz, 
as intended, but has the deleterious effect of increasing 
the radiation-pressure noise at lower frequencies. Conse- 
quently advanced detector sensitivity at almost all fre- 
quencies in the detection band is expected to be lim- 
ited by quantum noise. Qualitatively speaking, shot noise 
and radiation-pressure noise correspond to measurement 
noise and back action noise in quantum measurement 
theory — together they often impose the Standard Quan- 
tum Limit (SQL) to measurement accuracy . A correct 
modeling of the quantum noise of a GW interferometer 
should take into account correlations between the two 
types of noises, which may allow sub-SQL sensitivities to 
be achieved HSE|. 

The need for optical field simulation for gravitational- 
wave interferometer design has been addressed in the 
past with a variety of simulation tools, both in the fre- 
quency domain (e.g., twiddle and finesse 0) and 
in the time domain (e.g., the LIGO end-to-end simula- 
tion program Q). Although time-domain simulations can 
study issues associated with large mirror displacements 
and non-linear effects, e.g., the lock acquisition of the 
interferometer, they are computationally costly; in addi- 
tion, full time-domain simulations are also less straight- 
forward to quantize. In order to study the performance 
of gravitational-wave detectors, it suffices to stay in the 
linear regime near the operation point. For such a linear 
problem, frequency-domain simulations are dramatically 
simpler than time-domain ones; it is straightforward to 
obtain frequency-domain transfer functions, and there- 
fore noise spectra. In addition, since the system is lin- 
ear, the propagation of quantum Heisenberg operators 
are identical to those of classical field amplitudes, there- 
fore it suffices to build an essentially classical propagator. 

In low-power situations where radiation-pressure- 
induccd mirror motion is negligible and no non-linear op- 



tical elements (e.g., squeezers) are used, when linearizing 
over mirror displacements, propagation of electromag- 
netic fields at different frequencies are independent, and 
therefore the transfer functions can be established for 
each different frequency separately One only needs to 
take into account that, for the inputs to this linear sys- 
tem: (i) mirror motion (with frequency f2) creates phase 
modulation of the carrier, which is equivalent to generat- 
ing two equally spaced sidebands on the carrier frequency 
(at w±f2, where to is the carrier frequency and we denote 
lu + f2 and to — f2 as the upper and lower sidebands, re- 
spectively) with opposite amplitudes, and that (ii) laser 
noise can usually be decomposed into amplitude noise 
and phase noise, with the former contributing equally to 
the upper and lower sidebands, and the latter oppositely. 
These considerations have been the conceptual founda- 
tions of previous frequency-domain simulation programs. 

For high-power interferometers, the above strategy will 
have to be modified: the radiation-pressure forces acting 
on the mirrors, at frequency f2, depend on both upper 
and lower sideband fields; the induced mirror motion will 
again contribute to both sidebands — this makes it nec- 
essary to propagate pairs of upper and lower sidebands 
simultaneously. The mathematical formalism most con- 
venient for this problem, at least in the case of only 
one carrier frequency, is the Caves-Schumaker two-photon 
formalism |9lllCll|. In this paper, we adopt this formalism 
and present a mathematical framework for calculating 
the propagation of fields in an arbitrary optical system 
that includes the dynamical response of the mirrors to the 
light field. Namely, we divide complex interferometers 
into inter-connected elementary subsystems, and provide 
a general procedure for building a set of linear equations 
for all optical fields propagating between these systems 
- based on each individual system's input-output rela- 
tion, i.e., transformation matrices relating output fields 
to input ones and the incoming GW. We also describe 
the way in which these subsystems are connected to each 
other. Solving these equations will provide us with the 
optical fields, in terms of vacuum fluctuations entering 
the system from open ports, laser noise, and incom- 
ing GWs. While this mathematical framework, and the 
resulting numerical simulation tool, were developed to 
model quantum correlation effects in gravitational-wave 




FIG. 1: A sample configuration is shown. A beam block 
is connected to a mirror, which is in turn connected to a 
detector. Input fields incident on the mirror, a and d, are 
related to the output fields, b and c, by matrix operators 
derived in Sections ITT1 ITTT1 and ITvl 



interferometers, the method is general and can be used 
in any system where optical fields couple to mechanical 
oscillation modes. 

The paper is organized as follows: In Sec. [D] we intro- 
duce the mathematical framework for the simulation, and 
illustrate it with a simple example; in Sec. lIIII we provide 
input-output relations of basic optical elements that may 
be present in a laser interferometer, ignoring radiation- 
pressure effects and the presence of gravitational waves — 
by re-formatting well-known results in optics; in Sec. lIVI 
we take radiation-pressure-induced mirror motion into 
account, and provide input-output relations for movable 
mirrors and beamsplitters (up to linear order in mirror 
motion) , which have not been obtained before in the most 
general form; in Sec. we take into account the pres- 
ence of GWs by introducing modulation of cavity lengths, 
and treat the corresponding effect on light propagation 
up to linear order in L/Aqw (with L the length of the 
interferometer). In Sec. IVII the formulation is applied to 
a novel interferometer designed to extract squeezed vac- 
uum states that are created by a strong opto-mechanical 
coupling; and, finally, conclusions are summarized in Sec- 
tion En 



II. MATHEMATICAL FRAMEWORK 

A. General Prescription 

As mentioned above, the presence of opto-mechanical 
coupling dictates that we propagate the upper and lower 
sidebands simultaneously, which means that for each fre- 
quency f2, we will have to work with the two-dimensional 
linear space spanned by t he upper [a{oj + f2)] and lower 
[a(cj — f2)] sidebands [ly. Within the two-photon for- 
malism, developed by Schumaker and Caves 0, Il0l | , and 
outlined in Appendix lAl below, instead of a(uj ± Q), the 
two quadrature fields ai.2(fi) are chosen as the basis vec- 
tors. For simplicity of notation, we generally denote 

•-(2) <» 

and suppress the dependence of a on Q. 



We consider optomechanical systems formed by the fol- 
lowing elementary subsystems: movable mirrors, beam- 
splitters and free space propagators. We will also in- 
clude a "linear squeezer" , which turns an ordinary vac- 
uum state into a two-mode squeezed field with arbitrary 
squeeze factor and squeeze angle. Auxiliary to these 
optical elements, we introduce the beam block and the 
photodetector to deal with open ports which are cither 
left undetected and those detected with unit quantum 
efficiency; we also introduce the laser as an optical el- 
ement, which injects monochromatic carrier light and 
laser noise into the interferometer. Quadrature optical 
fields undergo linear transformations when propagating 
through such elementary systems, and quadrature fields 
with different fl's propagate independently from each 
other. These linear transformations are described math- 
ematically by the input-output relation, namely a set of 
equations relating the output fields to the input ones, in- 
cluding vacuum fluctuations, the carrier laser and laser 
fields, as well as to incoming GWs. We provide these 
input-output relations in Sees. IlIlllVl 

However, we note that propagation of sideband 
quadratures (Q ^ 0), although independent from each 
other, all depend on the propagation of the carrier 
quadratures (fl = 0), i.e., the amplitude and phase of 
the carrier incident on each subsystem. Fortunately, the 
propagation of the carrier is not affected by that of the 
sidebands, and can be carried out independently at the 
beginning. This said, we begin to formulate our general 
method of simulation. 

We build the following system of linear equations (for 
each sideband frequency Q) 
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where , i = l,...,N arc the N quadrature fields 
(each of them a two-dimensional vector) propagating in 
every part of the system, u^, i = 1, . . . N arc N gen- 
eralized input quadrature fields (each of them again a 
two-dimensional vector). The Mij, i,j = 1,...N are 
2x2 matrices which depend on the details of the optical 
system, and the can be written schematically as 

u W = v W+iW+hW/i > (3) 

where arises from vacuum fluctuations entering from 
the detection port or other lossy ports ( Sees . mTH 11111 and 
ITV|1. 1« from the laser (Sec. Ell, and H«/i from GW- 
induced phase modulation, with h the GW amplitude 
(Sec. [Vj; depending on the location of this generalized 
input field, some or all of the above three contributions 
could also be zero. Henceforth in the paper, we shall 
consider each pair of quadrature fields as one object. In- 
verting the matrix Mij will give a^ in terms of vS l \ and 
hence all of the necessary transfer functions. 

Now let us provide a universal prescription for con- 
structing Eq. J2J), suitable for modelling generic systems. 
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Wc break this procedure into two steps: 



To comply with the format of Eq. , we write 



1. Suppose we have n elementary subsystems men- 
tioned above, with the fc th subsystem having p^ 
ports. The entire system will then have P = 
J2k=iPk ports. Because we formally include beam 
blocks and photodctcctors as subsystems, none of 
our ports will be formally open, i.e., left uncon- 
nected to some other port. This means that we 
have P/2 pairs of connections. For each pair of 
connections, we have two fields, one propagating 
in each direction. This means we have a total of 
P fields (each in turn has two quadrature compo- 
nents). 

2. For each system k, with p^. ports, we also have pk 
input fields and pk output fields, and therefore the 
input-output relation will provide us pk equations. 
All subsystems together will then provide us with 
P equations (each with two components), exactly 
the number needed. 



B. Example with the input-output relation of 
beam blocks, photodetectors and lasers 

Next we illustrate the generic construction procedure 
with a simple example, which also clarifies the formal 
roles of beam blocks, photodetectors, and lasers. We 
first propagate fields between three basic elements of an 
optical train: a beam block, a partially reflecting mirror, 
and a photodetector. Referring to Fig.^ the beam block 
is connected to the mirror, which is in turn connected to 
a detector. For simplicity, we assume that the mirror is 
lossless and fixed in position. 

As a first step, we identify the fields in consideration. 
The beam block and the photodetector are 1-port sys- 
tems, the mirror is a 2-port system; we have a total of 
4 ports, and 4/2 = 2 connections. There are two fields 
associated with each connection; we label them a, b, and 
c, d, respectively, as done in Fig.^ Note that each field 
in turn has two quadrature components, so the system is 
8-dimensional, and we need 8 scalar equations. 

Now we have to provide the input-output relations for 
each object. For the mirror with amplitude reflectivity 
p and transmissivity r, and neglecting radiation pressure 
effects, we have 




(4) 



Note that Eq. (@J contains 4 scalar equations, and that 
p and r are really 2x2 scalar matricies, pi, and rl (this 
is true because our mirror does not mix quadratures) — 
we have suppressed the identity matrix I for simplicity. 
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For the beam block and the photodetector, they really 
are placeholders for physically open ports. Their input- 
output relation is simply that the output fields from 
them are vacuum fluctuations (independent from the in- 
put fields): 
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Here we assume implicitly that the photodetector is de- 
tecting the field c with unit quantum efficiency. In order 
to model imperfect photodetectors, we could add a mir- 
ror with zero reflectivity and non-zero loss in front of the 
ideal photodetector. 

Combining Eqs. JSJ) and ©, we have 



f -1 

-p -1 

T — 1 

V o o o 



\ 



I J 



b 

c 



( -v (1) \ 





(7) 



M 



which are the 8 scalar equations we need. Inverting M. 
will give us each of the propagating fields in terms of the 
input vacuum fields. 

Now suppose the beam block is replaced by a laser 
source, coupled to the spatial mode of a field, then we 
only need to replace the vacuum field in Eqs. © and 
by the laser field, 1 (1) : at Cl = 0, carrier quadratures, 
while at CI ^ 0, it gives the laser noises. 

Here we note that all diagonal elements of A4 are equal 
to —1 — this is in fact not a coincidence, but a universal 
feature of our construction procedure. In order to un- 
derstand this, we need to realize that every field a( fe ) is 
the output field of exactly one subsystem. In the input- 
output relation of that unique subsystem, there is exactly 
one line that relates a^ fc ) to the input fields of this sub- 
system, which reads: 
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This equation corresponds to, after moving aS k ^ to the 
right-hand side of the equation, moving any non-a 1 
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1, . . . , N terms to the left-hand side, and swapping left 
and right, 
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It is obvious that the lines of equation found by this way 
for different a^'s will be different. As a consequence, 
we can arrange to have the line corresponding to a' fc ' 
appear on the k th row of Ai, and thus have all its diagonal 
elements equal to —I. 




Free space 




III. MATRICES FOR STATIC OPTICAL 
ELEMENTS 

In this section, we derive the matrices for some stan- 
dard objects used in simulating quantum noise in a 
gravitational-wave interferometer. Here we neglect ra- 
diation pressure effects and the presence of gravitational 
waves (they will be dealt with in Sees. II VI andlVl respec- 
tively). As a consequence, our derivation only involves 
some re- formatting of previously well-known results. 



A. Mirrors 

Field transformations due to a mirror were introduced 
in the example of Section [Q] The transformation ma- 
trix for a lossless mirror is given in Eq. (@J. We now 
derive more complete equations for the mirror that in- 
clude losses. We ascribe a power loss A to the mirror 
in Fig. ^such that p 2 + t 2 + A = 1. The introduction 
of losses gives rise to an additional vacuum field of am- 
plitude J A/(l — A) that is added to each input of the 
mirror. The (I — A) 1 factor accounts for part of the 
loss field being lost to the mirror. This can be verified 
by having shot-noise-limited fields, a and d, incident on 
the mirror. The field returning to the beam block 
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must also be at the shot noise level, such that 
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The new equations governing the mirror are 
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where v*- 3 ) and v^ 4 - 1 are the vacuum fluctuations that en- 
ter due to the presence of loss. 

Equation (|12|l may be rewritten as 
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FIG. 2: The fields entering and exiting a region of free space 
are shown. Propagation operators are characterized by the 
propagation distance (and orientation relative to the source 
polarization, in the case of the GW signal). 



where 



43)/ 



v (4)' = 



1-A 



1-A 



TV^+pv' 4 ' 



(14) 
(15) 



v^ 3 )' and v^ 4 )' are uncorrelated vacuum fields in this rep- 
resentation. We can subsequently write the mirror's con- 
tribution to Eq. J5J as 
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This method may also be used to inject losses in beam- 
splitters or cavities. 



B. Free space propagation 

Since optical cavities are present in virtually all optical 
configurations of gravitational- wave interferometers, we 
must give a transformation matrix for them as an element 
of our arbitrary optical train. To do so we introduce an 
operator to transform the field as it propagates through 
free space between any two other optical elements (in 
the case of an optical cavity, these would be mirrors). 
Using the convention of Fig. [21 the matrix for propaga- 
tion through a length L transforms input fields a and d 
according to 



where the matrix for the propagator is 
M Prop = e** 

\K e , 

Here 
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FIG. 3: Treating the beamsplitter as a four-port device, def- 
initions for the fields, including sign conventions, are shown. 



D. Correlators 

The correlator module of the simulation allows for the 
inclusion of squeezed light or vacuum fields in the inter- 
ferometer. It is essentially a one-way device: only fields 
entering from one direction are transformed; fields enter- 
ing from the other direction pass through the correlator 
unmodified. Taking a to be the input field, the field at 
the output of the correlator, b, is defined by 

b = S(r,0)a, (24) 

where S (r, 6) is the squeeze operator with squeeze factor 
r and squeeze angle 6: 

S(r,6) = 

(cosh r + sinh r cos 26 sinh r sin 26 \ 

■ (25) 
sinh r sin 26 cosh r — sinh r cos 26 I 



are the one-way phase shift on the carrier light at fre- 
quency, to, and on modulation sidebands at frequency, O, 
respectively, and 

(cos — sin \ 
(21) 
sin cos / 

is the rotation operator on quadrature fields. 

C. Beamsplitters 

Another essential optical element of an interferome- 
ter is the beamsplitter. We consider a beamsplitter with 
amplitude reflectivity and transmissivity, p and r, re- 
spectively. The beamsplitter transforms the input fields, 
shown in Figure |3 according to the matrix equation 
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where 
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M BS = . (23) 
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In presence of optical loss, assuming p 2 + t 2 + A = 1, 
and going through similar arguments to Sec. IIII Al we 
simply add a column vector of vacuum fields — y/Av^ 
(i = 1, 2, 3, 4) onto the right-hand side of Eq. JUJ). 



IV. RADIATION PRESSURE 

Radiation pressure plays an important role in inter- 
ferometers operating close to or beyond the SQL, since 
quantum back-action noise must be taken into account. 
Moreover, radiation-pressure effects can also modify the 
dynamics of these interferometers ||| . Sideband quadra- 
ture fields create amplitude modulations to the carrier 
field, and the associated power modulation drives the mo- 
tion of optical elements, which, in turn, phase modulates 
the carrier, thereby creating sideband quadrature fields. 

Details of this sideband-to-sideband conversion depend 
on the phases (this determines which quadrature gets 
converted into which) and amplitudes (this determines 
the conversion strength) of the carrier field propagating 
in different parts of the interferometer. Therefore, it is 
necessary to separate the fields into carrier (f2 = 0) and 
sideband (f2 ^ 0) components at this point. The radia- 
tion pressure force due to the carrier field itself is a time 
independent force and can be ignored (in reality they will 
be balanced by a static force exerted on the optical ele- 
ments, e.g., the pendulum restoring force on a suspended 
mirror). The effect of interest is the time-dependent part 
of the force, due to sideband components, which will be 
the subject of this section. As a foundation, we must 
first of all calculate the phase and amplitude of the car- 
rier fields at each location. But this we can already do 
by building the general equation J2J out of input-output 
relations of static optical elements, which have already 
been derived in Sec. IIIII and solving it. 

Before incorporating radiation pressure into the treat- 
ment of specific systems, let us study the electromag- 
netic momentum flux carried by optical fields in the two- 
photon formalism. In quadrature representation, we de- 
compose the total quadrature field Ej otal (here Ej can 
be a, b, c or d for the configuration in Fig. ^) into the 
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following two terms: 
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(26) 



The monochromatic carrier field in Eq. 1)26(1 can be 
written more explicitly in terms of power Ij, phase 9j 
and effective beam area A as 
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while the sideband field can be written as an integral over 
all sideband frequencies: 
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The total momentum flow carried by the field is 
A 
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Removing the static (dc) and optical frequency (u) com- 
ponents, the Fourier transform of the timc-avcraged (over 
a time scale much shorter than the GW period, but much 
longer than 1/uj) ac momentum flow carried by this field 
is 



i>-(n) 



where we have defined 
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as the carrier quadrature field, and j(f2) is the sideband 
component at angular frequency £1. 

In the remainder of this Section we derive explicit 
input-output relations for mirrors and beamsplitters, in- 
cluding radiation pressure effects. Our results will be 
more general than previously obtained results by allow- 
ing the carrier fields incident from different ports to have 
different phases. 



A. Mirrors 

Let us once again consider the mirror in Fig. As- 
suming that the mirror behaves as a free particle with 
mass M when no radiation-pressure forces are exerted 
(valid for suspended mirrors when frequencies greater 
than the pendulum resonant frequency are considered), 
the Fourier transform for the equation of motion for the 
mirror is 
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where X is the displacement of the mirror induced by all 
the sideband fields (X is positive to the left in Fig. ^ and 
the j refer to a, b, c,d). The summation is performed 
over all the fields entering and exiting the mirror; the 
coefficients rj a = rjb = — 1 and ?/ c = r]d = 1 account for 
the directions of propagation. The displacement of the 
mirror due to the radiation pressure forces, X, can be 
written explicitly as [see Eq. 13UI) ] 
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Given a (time-dependent) displacement X{t) of the mir- 
ror, the input-output relation can be written as (if X <C 
c) 
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c in the argument of iJj otal for the j-th field is the speed 
of light and should be distinguished from c in the sub- 
script of Ej° tai , which refers to the field c. In quadrature 
representation, to leading order in X and in the sideband 
field amplitudes, we have, 
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, for v = 













Here * refers to a rotation by 7r/2, as described by 
— Re=ir/2 m Eq. I|35|l . Accordingly, for any quadrature 
field v, we define 



(36) 



Equation (|35() implies that time delays, or phase mod- 
ulations, create sideband quadratures orthogonal to the 
carrier, as illustrated in terms of phasors in Fig. 21 The 
sideband part, i.e., the ac components in Eqs. I|34a|) 
and l|34bj) . can be obtained using Eq. (|3*5|) : 
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Inserting Eq. (J33J) into Eq. lj37|l gives 
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FIG. 4: Here we show that the phase modulation sideband 
generated by the radiation pressure force is perpendicular to 
the carrier field, which is why the generated signal has a D* 
dependence. 



where 



IT 



2 p u> 

m n 2 c 2 



(39) 



is a quantity with units of inverse power or W _1 . [For 
lossy mirrors with p 2 + t 2 + A = 1, we simply insert 
a column vector — \/Av^ , i = 1,2 onto the right-hand 
sides of Eq. O and ®, Cf. Sec. lrlLAl ] 

To solve for b and c, the matrix on the left hand side 
of Eq. I)38|) must be inverted. It is straightforward to find 
a complete set of eigenvectors for this matrix, they are: 



Physically, this resonance comes about because the side- 
band fields generated by mirror motion can exert radia- 
tion pressure back onto the mirror. Let us for a moment 
consider classical motion of the mirror. As was men- 
tioned after Eq. (|37|l . for any given input carrier field, the 
sideband field generated upon reflection from the moving 
mirror is it/2 phase shifted relative to the input carrier, 
so the sideband will not beat with the reflected carrier to 
induce any force on the mirror [see Eq. (|30|l ] — force can 
only be induced by beating this motion-induced sideband 
field with the transmitted carrier, which must have non- 
zero amplitude and must have a phase difference other 
than 7r/2 relative to the sideband. This explains why the 
resonant frequency vanishes if either p = or r = (no 
reflected or transmitted field), or if 6 a — Od = Nir (no 
phase difference between the two input fields). 

When the two input carrier fields, D a and D^, have the 
same phase (or differ by Nir) , the phasors corresponding 
to D Q , Db, D c and D<j all become parallel to each other. 
This is true for almost all interferometers that have been 
treated explicitly analytically. This case is rather special 
from a mathematical point of view, since the matrix we 
are inverting does not have a complete set of eigenvectors. 
Fortunately, the inverse is just 
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Since the first three vectors are orthogonal to 
( — D J ) , the three corresponding eigenvalues are 
Ai = A2 = A3 = 1; the last eigenvalue is 
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Inverting the eigenvalue A4 yields a pair of resonant fre- 
quencies at 
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(This identity originates from the fact that the sideband 
field is orthogonal to the carrier field about which it is 
generated.) Using this fact, we can further simplify the 
input-output relation to 
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[Here for simplicity we have assumed the mirror to be lossless.] In practice, although Eq. 1)38(1 does not give the 
output fields b and c explicitly in terms of the input fields a and d, it can be incorporated to the matrix M. (and into 
, in presence of optical losses) without any trouble [cf. Eq. |J2J]: its inversion will take place automatically when 
M^ 1 is calculated. [However, doing so will make it impossible to have —1 all along the diagonal of M.] Alternatively, 
the variable X may be added to our system of variables, with Eq. Q.'S3[1 providing the additional equation necessary. 
The equations governing a mirror may then be replaced with Eq. I|37|l to include the dependence on X. In this way, 
the —1 diagonal components are preserved, without the need to invert additional matrices. 
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B. Beamsplitter 



Referring to the fields shown in Fig. |3| the displacement due to radiation pressure forces on a beamsplitter (normal 
to its reflective face) is 



X X + Xy 



1 



Am 2 V 2c 2 



T>T, D 



D„ D 





(46) 



where X x is the displacement along the x-axis and X y is the displacement along the y-axis. Similar to the case of a 
cavity mirror, this motion induces phase fluctuations on the impinging fields upon reflection, and introduces additional 
terms in the input-output relation. Following a procedure similar to the one with which we obtain Eq. I|37|) , we get 











= M bs 




(l 











cy/hui 



Inserting Eq. (|46|l into Eq. I|47|l gives 



D* 

V D // 



D: D 



D 













= 















M BS - ? 



D* 

V D // 



D* 

V D // 



(47) 



(48) 



Equation (|48|l is quite similar in nature to Eq. 1)38(1: optical losses can also be incorporated in a similar fashion, by 
adding — VAvW, i = 1,2,3,4 on to its right-hand side, where p 2 + r 2 + A = 1. Again, in the generic case where 



D 



^0, 



(49) 



the matrix on the LHS of Eq. I|48|l has eight linearly independent eigenvectors, of which seven have unit eigenvalue, 
while the eighth has 



= l + rn(DTDS + D^D fe *) 



V D /7 

eigem 



D 6 

V D /7 



= i 



4/5TW0 



Vlfldsm(0f 



Mfl 2 c 2 

which corresponds to an opto-mechanical resonance at angular frequency 

Apruj 



b smiVh - V b/ 



±n BS = ±t - 



Mc 2 



y/hT b sm(6 h - e b ) + y/IJhsm(e f - e d ) 



1/2 



In the special case of 



Df D 



D 



D* 
D* 

V D // 



(50) 



(51) 



(52) 



i.e., all input carrier fields are in phase with each other (modulo 7r) wc get 



e 

Vg 









D* 


M BS - pH 









( DJ T>1 D 





if D h II D rf II D ( 



(53) 



For simplicity, we assume the beamsplitter to be loss- 
less in the above equation. This is particularly true for 
the beamsplitter in Michelson- and Sagnac-type GW in- 
terferometers Similar to the case of the mirror, for 
the purposes of simulation, we incorporate the position 
of the beamsplitter as an additional variable in M, in 
order to preserve the —1 diagonal elements and to avoid 
the inversion of additional matrices. 



three components: the gravitational-wave signal, classi- 
cal laser noise, and noise due to vacuum fluctuations in 
the detected mode. The outgoing field being detected, b, 
has the general form [see Eqs. (J2J) and (pj]: 



i 



(56) 



V. GRAVITATIONAL WAVE SIGNAL AND 
THE OUTPUT FIELD 

A. GW contribution 

In our set of optical elements, only optical cavities have 
significant propagation distances, so we model the effect 
of GWs by introducing a phase shift to the carrier light as 
it passes through a cavity. To calculate the propagation 
of these fields, all that must be done is to add a source 
term in the equation governing the cavity. Refering to 
the fields in Fig. [21 the cavity field becomes 



c = 



e *Rea-»/-^D; 



Re 




(54) 



where h is the Fourier transform of the GW amplitude. 
An /i-dependent term is also added to the equation relat- 
ing b and d using Ti* d in place of D*. The parameter r\ 
takes values from — 1 to 1 depending on the orientation 
of the cavity and the polarization state of the incoming 
GW. For example, for a linearly polarized incoming GW, 
and for an optimally aligned Michelson interferometer, 
we have r\ = 1 for one and —1 for the other. 

It is straightforward to incorporate Eq. I)54|l into the 
general equation Eq. Q. In particular, the term contain- 
ing h on RHS contributes to the GW part of the general 
input field u, i.e., to the third term of Eq. ©, with 



H = -7] 



uiL 



2c\/Ti<jj 



D* 



(55) 



B. Photodetection: signal and noise 

For our purposes, the photodetector serves two roles: 
first, it represents an open port, from which vacuum fluc- 
tuations enter the interferometer; second, it determines 
the measurement point. For the former, the input-output 
relation of a photodetector, as it contributes to the ma- 
trix M and the generalized input vector u«, is trivial 
and has been discussed in Sec. Ill Bl Here we focus on 
the latter. At zero frequency, there is only contribution 
to b from the carrier laser, while at non-zero sideband fre- 
quencies, the detected fields at a photodetector comprise 



The summation is performed over all fields. We note that 
contributions to vW exist only for fields that emerge from 
beam blocks or lossy optical elements, those to l^-* exist 
only for the field that emerges from the laser, and those 
to J3y> only for fields that emerge from cavities. 

We suppose homodyne detection at quadrature angle, 
is performed such that the measured field is 



b\ cos £ + 62 sin (. 



(57) 



For a complete simulation, C should be the phase of the 
carrier that emerges at this port. However, in theoretical 
studies, we could also assign another value to C, assuming 
that the local-oscillator phase is modified by some other 
means the simulation does not address. 

For the detected field, the quantum noise spectral den- 
sity is (see, e.g., Sec. Ill of Ref. Q) 



[cosC sin (]T M S Vl T b \ 



cos( 
sin £ 



(58) 



Because is always proportional to a vacuum field, we 
have used S Vi to denote the noise spectral density which 
is identical for all its quadratures. Here we have added 
the power of different loss contributions, since we as- 
sume the vacuum fields to be independent to each other. 
In general, laser noise is neither quantum-limited, nor 
are the magnitudes of phase and amplitude fluctuations 
equal; there could also be correlations between the laser 
amplitude and phase noise, even as the laser field enters 
the system. Taking these into account, we have a laser 
noise spectral density of 



(Nl) b = [cosC sin( ]T bl S L T b \ 



cos£ 
sinC 



where I corresponds to the input laser field, and 



Si 2 



S12 
S22 



(59) 



(60) 



describes noise of the laser as it first enters the interfer- 
ometer, with jSii(22) being the noise spectral density of 
the first (second) quadrature, and S12 the cross spectral 
density between the two quadratures. [In the usual con- 
vention of having input laser in the first quadrature, 1 
corresponds to the amplitude quadrature, hence ampli- 
tude noise, while 2 corresponds to the phase quadrature, 
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hence phase noise.] The transfer function for the GW 
signal is 



H b = J2[ C0S C sinC]T bt H 



(i) 



(61) 



Note that GW contributions from different parts of the 
system add up coherently. The displacement (strain) 
noise spectral density from quantum noise is then given 
by 



Sh 



N 2 Q + Nj 



(62) 



VI. 



APPLICATION TO A COMPLEX 
INTERFEROMETER 



The mathematical formulation described in Sections ITTI 
through Ivl was encoded into a simulation program writ- 
ten in C++. In this section we describe tests of the 
simulation code for a complex interferometer configura- 
tion, where the simulation results were compared with 
analytic calculations. 

The interferometer configuration is shown in Fig. [5J 
and in Fig. we show fields propagating in the interfer- 
ometer as well as modes of motion of the mirrors. The 
interferometer is similar to that used in GW detection: 
a Michclson interferometer with Fabry-Perot cavities in 
each arm. All the mirrors of the interferometer are sus- 
pended as pendulums. Power-recycling |T3 | is optional 
and is not included here. The configuration shown has 
a few unusual features compared with a conventional in- 
terferometer, however. First, the end mirrors of the arm 
cavities are a common suspended object, coated with a 
high-reflectivity coating on both surfaces and assumed to 
have an opaque substrate. Second, this cavity end mir- 
ror object is very light, with a typical mass of 1 g, and 
is suspended as a pendulum with resonant frequency of 

fold mirror 

^cavity input mirror (250 gm) 

_ shared end mirror (1 gm) 
cavity input mirror 

fold mirror 




about 1 Hz. All remaining optics are assumed to have 
a mass of 250 g, and are also suspended as pendulums 
with a resonant frequency of 1 Hz. Third, the cavities 
are detuned from resonance. 

Testing the simulation with this somewhat unconven- 
tional interferometer configuration served two purposes: 
(i) It is the baseline design for an experiment to gener- 
ate squeezed states of the electromagnetic field, produced 
with radiation-pressure-induced optical forces in an in- 
terferometer with low-mass mirror oscillators and high 
stored power [Tflj : and (ii) the shared end mirror gives 
rise to unexpected dynamical effects that prove interest- 
ing and instructive to explore, and are relevant to other 
high-power interferometers, such as Advanced LIGO Q. 
We note that the shared end mirror has advantages in 
terms of mechanical stability and control system design, 
but the desired radiation-pressure effects can be realized 
by a configuration with two independent end mirrors as 
well. 



Parameter 



Symbol Value Units 



Light wavelength Ao 

End mirror mass m 

Input mirror mass M 
Input mirror transmission Ti 

Arm cavity finesse T 
Loss per bounce 

Arm cavity detuning (f> 

Input power Io 



1064 
1 

0.25 

4 x 10~ 4 
1.6 x 10 4 

5 x 10~ 6 
10~ 5 

1 



kg 



W 



BS reflectivity asymmetry A B s 0.01 
Michelson phase imbalance Aa M 
Michelson loss imbalance Ae M 
Input mirror mismatch At 5 x 10 -6 - 
Detuning mismatch A^ 10~ 7 Ao 

Arm cavity loss mismatch A e 2 x 10 -6 - 

TABLE I: Select interferometer parameters and their nominal 
values. 



e bandwidth (T< + T e )c/(4L) 

el bandwidth due to loss T e c/(4L) 

— A resonant frequency <f>c/L 

a characteristic quadrature rotation angle arctan(A/e) 

TABLE II: Quantities associated with the detuned arm cavi- 
ties . 



FIG. 5: Schematic of a an interferometer designed to extract 
ponderomotively squeezed light due to radiation-pressure- 
induced motion of the ultra-light shared mirror. Light from 
a highly intensity- and frequency-stabilized laser source is in- 
cident on the beamsplitter. High-finesse Fabry-Perot cavities 
in the arms of the Michelson interferometer are used to build 
up the carrier field incident on the end mirrors of the cavity, 
which are a single mechanical object. 



A. Ideal optical springs 

In this section we study analytically a crucial com- 
ponent of the interferometer design: the optical spring 
effect, especially in the case of two identical detuned 
cavities with a common end mirror. The input-output 
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FIG. 6: Optical fields propagating in the interferometer, and 
modes of motion of the mirrors. In particular, ckma.b and 
ema.b are artificial detunings and losses one can add to the 
two arms of the Michelson interferometer, respectively, see 
Sec. lVllfl for their significance. 

relation of this system can be obtained by carrying out 
our generic procedure analytically. In doing so, we ex- 
tend previous results in Refs. [a, [3 to include two new 
features. First, we consider motions of all three mirrors, 
with mass of the input mirrors different from that of the 
common end mirror. Second, in our system the carrier 
phases incident on mirrors are different; under such a cir- 



cumstance, formulas developed in Sec. II VI are non-trivial 
extensions to existing ones. 

In order to make results intuitively understandable, we 
consider only the ideal system, with the two input mir- 
rors completely identical, the common mirror perfectly 
reflective on both sides, the two cavities having exactly 
the same lengths, the carrier incident on both input mir- 
rors having equal amplitude and phase, and with a per- 
fect beamsplitter. We also ignore the free pendulum fre- 
quency, and consider the test masses to be free. Similar 
to previous studies, we assume a high-fincsse cavity and 
ignore the interaction between the motion of the input 
mirror and the carrier light outside the cavity. We retain 
terms only to the leading order in eL/c, XL/c and SlL/c, 
where L is the cavity length, c is the speed of light, SI is 
the sideband frequency, and (—A — ie) is the complex op- 
tical resonant frequency of the cavity with fixed mirrors 
[—A denotes the resonant frequency and e the bandwidth, 
defined in Tabic ITT1 and wc ignore end-mirror loss]. 



1. Differential Mode 

With the above assumptions, the differential optical 
mode couples only to differential modes of mirror motion: 
those with the two input masses moving such that xa = 
—xb = x-d, and arbitrary x m [see Fig. 0; such modes 
form a two-dimensional subspace of all possible motions 
of the three mirrors. In the ideal case, we only need to 
study this mode. The differential input-output relation 
is given by 




with 



C D = 



Rq, 



C d R_q, ( 1 ) 


+ s D 


V a2 / 





-(SI 2 - A 2 + e 2 )Sl 2 - Xl d 
-2 £ AS1 2 + 2et D 



,(o) , J") 



2eAQ 2 

-(S! 2 -A 2 + e 2 )Sl 2 -Ai. D 




Lh SQL 



(63) 



(64) 



and 



M D = £l 2 [(SI -He) 2 - A 2 ] + At D 



(65) 



Here i™ is the motion of a free end mirror with the 



same mass, if is the free differential motion of the input 
mirrors (x^ = —Xg = 4°'); ol = arctan(A/e) is the 
carrier phase at the end mirror. The carrier incident 
on the input mirrors has phase 0, the carrier inside the 
cavity, leaving the input mirror has phase a — (f>, while 
the carrier inside the cavity entering the input mirror 
has phase a + <fi. The quantity /igq L is the free-mass 



Standard Quantum Limit associated with the differential 
mode, given by 



= 2mM/(m + 2M) . (66) 



The quantity t D , defined by 

t D = t ' ( 67 ) 

measures the strength of optomechanical coupling [notice 
the dependence on carrier intensity I c and the inverse de- 
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pendence on the effective mass of the differential mode 
mechanical oscillator /x D ]. Roots of M D are the (com- 
plex) resonant frequencies of the coupled optomechanical 
system. From t D we define a characteristic frequency, 



A 2 ) 



(68) 



For systems with D <C e, the two resonances are well 
separated, and are given approximately by ±0 D [me- 
chanical frequency due to optical spring] and (±A — ie) 
[optical resonant frequency] , respectively — this is indeed 
the regime in which we construct our experiment. 

The differential optical mode couples to a two- 
dimensional subspace of all possible motions of the three 
mirrors. It is instructive to look at the motion of sepa- 
rate mirrors, in the regime of Q <C e, i.e. for sideband 
frequencies f2 well within the linewidth of the cavities: 



J n 



e 2 -n 2 
e 2 



A 2 



1 

e 2 



-n 2 



A 2 e? 



A 2 - 

A 2 e 2 

'A 2 + l A 2 + l 



Q 2 



Jo) 

JO) 



(69) 



Here we have defined A 2 = 2M/m. From Eq. (|69|) . we 
conclude immediately that 



n 2 



e 2 ~n 2 



a:< > + x {0) 



(70) 



This change in response is exactly what happens when a 
free test particle is connected to a spring with mechanical 
resonant frequency D . Equation (|70|l reveals a crucial 
advantage of the optical spring — that the response of the 
cavity length to external disturbances (driven by seismic 
and/or thermal forces, e.g.) is greatly suppressed from 
the corresponding value for free-mass systems. Theoreti- 
cally, this suppression is present even when a mechanical 
spring is used. However, mechanical springs introduce 
thermal noise, usually of much higher magnitude due to 
the intrinsic mechanical loss [l4l ll5j] . 

It is interesting to notice that the suppression of total 
cavity length fluctuations is achieved collectively by the 
end mirror and the input mirror. As we see from Eq. 1)69(1 . 
[in the case of large A] , the motion of the end mirrors x m 
is suppressed from its free mass value by the factor in 
Eq. H70|) , while the motion of the input mirrors x D is not 
influenced by the spring, since it is relatively massive. 
Fortunately, through the (1,2) component of the matrix 
on the RHS of Eq. 1|69[) . this motion of the input mirror 
is imposed onto the end mirror with opposite sign, again 
suppressing the total cavity length fluctuations. 

Now let us restrict ourselves to the regime of il < D < 
e, and study the quantum fluctuations and classical com- 
ponent of the output field (due to classical disturbances 
to the mirrors). As we shall sec shortly, this regime has 
two crucial features: (i) the response of the output field 



to im + iEd , and thus length fluctuations due to seismic 
and thermal noise, are greatly suppressed by the optical 
spring and (ii) the output squeezed state is frequency- 
independent. 

For quantum fluctuations, we have 



C D 



-1 
2e/A -1 



(71) 



which is frequency-independent. It is straightforward to 
derive that the quantum noise spectrum in the = 
bi cos C + &2 sin £ quadrature [Cf. Eq. I|58|l ] : 



1 



A 2 



- 2 



+ ^ cos(2C - 3a) 



(72) 



In particular, terms in e/A are associated with squeezing, 
where the constant power squeeze factor e 2q (q > 0) is 
given by 



(73) 



-2q\ 



sinhg = |e/A| . 

The minimum noise spectral density (S^ = 
reached at £ = 3a/2, while at £ = a and 2a the noise 
spectrum is equal to the vacuum level (S^ = 1). Values 
of e/A corresponding to several power squeeze factors are 
listed in Table ITTll As shown, e and A will not differ by 
a factor of more than ~ 2, for typically desired squeeze 
factors. 

Now for the classical component, given by the second 
term in Eq. (|63|) . we have 



1 



-Rn 



Lh SQL 



f£2=_e 

e 2 a 



sin 2a 
- cos 2a 



(74) 



This means the entire signal due to differential displace- 
ment if + x^ is in the single quadrature £ = 2a + ft/ 2, 

and there is no x„ + x^ signal in the £ — 2a quadra- 
ture. Interestingly, the quantum noise in this quadrature 
is right at vacuum level. In addition, since /ig QL oc 1/fi, 



the response of b^ to proportional to Q~ at 

this regime - therefore not only the motion, but also the 
output field, has a suppressed response to thermal and 
seismic noises. Note here that the suppression factor is 
proportional to \JT~ C (since 9 U oc ^JTZ oc \/T^) - because 
motion is suppressed by I c , while the optical sensing of 
mirror motion is enhanced by J~T C . Now suppose we in- 
troduce a noisy force which induces a spectral density 
on a free mass, then the output classical noise will 

be 



(o) 



S, 



4^Ysin 2 (C-2a) 



(75) 



At the minimum quantum noise quadrature, £ = 3a/2, 
we have 



D 3a/2 



2e 
: 0.6 



1 



n 2 



Va 2 ^ 

5 ;r 



o 2 



e 2 D L 2 (h- v 



© 2 d L 2 (h- 



(76) 
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where the inequality is obtained by taking maximum over 
all e and A. We note that because of the suppression 
factor f2 2 /0 2 , the classical noise S% can be much higher 
than the free-mass Standard Quantum Limit while still 
allowing the interferometer to generate squeezed vacuum! 



Squeeze Factor (dB) 


3 7 10 20 


e/X 


0.58 1.13 1.42 2.12 



TABLE III: Relationship between power squeeze factor and 
e/X, see Eq. lf?5jl . 



xa = xb = x c . This mode is irrelevant to an ideal 
interferometer with identical arms and perfect contrast. 
In reality, however, the common mode will influence the 
output via couplings induced by differences (mismatch) 
between the two cavities, for example. Such effects can 
be quite important near the common-mode optomechan- 
ical resonance. 

The input-output relation of the common mode, simi- 
lar to that of the differential mode [cf . Eq. (|63Jl ] , is given 
by: 



2. Common Mode 




CcR— a 




S a X 



(0) 



(77) 



We now consider the common optical mode, which cou- 
ples with motion of the input mirrors corresponding to 



with [cf. Eq. 



J 



-(O 2 



-A 2 +e 2 )f2 2 
-2eAf2 2 + 2et c 



Air 



(n 2 - x 2 + e 2 )n 2 - \l 



r 




(78) 



and [cf. Eq. (gSJ)] 

M c = n 2 [(n 



+ ief - X 2 ] +Xl c 



(79) 



/igq L , the SQL associated with the common mode, is 
given by [cf. Eq. i|fJrJ|> ] 



"■SQL 



2h 



fi c = 2M . 



(80) 



n D n 2 L 2 ' 

The quantity i c is given by [cf. Eq. (|rJ7|) 

8loqI c 
ji c Lc 

For the common mode, we have a optomechanical reso- 
nant frequency of [cf. Eq. I|68|l ] 



(81) 



e, 



N / ic A/(e 2 + A 2 ) , if9 c « 



(82) 




This frequency is in general much lower than its 
differential-mode counterpart, with 



(83) 



B. Laser coupling to the antisymmetric port due to 
mismatch 

Mismatch between the optical parameters of the two 
arm cavities, as well as imbalance in the beamsplitter re- 
flection/transmission ratio and imperfect contrast of the 



6 

AA 

Acj-m 
A BS 
Ae M 



eX 



e 2 + A 2 
e 

Ve 2 + X 2 
eX 

e 2 + A 2 



2a + ty/2 
a 

2a + tt/2 
2a + 7r/2 

2q 



TABLE IV: Transfer function from carrier to differential 
output [see Eq. (I91t ]. in the leading-order approximation. 
The same coefficients apply to phase-noise coupling, i.e., 
iV^ = C(fc), <p^ — <pfk)i m the low- frequency regime [see 
Eq. JSSJ]. 



Michelson interferometer, can couple the carrier and also 
the noise sidebands on the laser to the differential detec- 
tion port. For each arm, A and B, we denote the true 
value of the kth quantity by its nominal value plus con- 
tributions due to imperfections, i.e. 



X 



(k)A, 



= x (k) ± 1a; 



(fe) • 



(84) 



Here the index k refers to the type of imperfection being 
considered. The beamsplitter asymmetry is character- 
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Ae 

e 


e 2[ e 2 (n 2 +e 2 )2+4A 2 e 4]l/2 £ 3 2A q2 

A(e 2 + A2)(n2_ e 2) A(e 2 +A 2 ) " a aiCtan e(n2+e2) 


Ael 

e 


e 2 e 2 
AVe 2 + A 2 AVe 2 + A 2 


AA 


e [(A 2 fi 2 -e 2 e 2 ) 2 +46 2 A 2 e^ 1/2 £ A „ 2 e A6 2 

A(e 2 + A 2 )(fi 2 - 6 2 ) e 2 + A 2 " a ' arCtan A 2 ^ 2 - e 2 9 2 


Aq m 


[A 2 (« 2 -e 2 ) 2 + 4 e 2 e 4 c ] 1/2 1 o t 2ee 2 
A(n 2 -e 2 ) 1 2q arctan A(n 2 -e 2 ) 


A bs 


2e« 2 2e 

A(n 2 - e 2 ) a 


Ae M 




~[(e 2 + A 2 )fi 2 - (2e 2 + A 2 )9 2 ] 2 + e 2 A 2 9 4 " 


1/2 

Ve 2 + A 2 eX& 2 c 


2AVe 2 + A 2 (« 2 -e 2 ) 2A " (e 2 + A 2 )fi 2 - (2e 2 + A 2 )9 2 



TABLE V: Laser amplitude noise coupling into the dark port, in the leading-order approximation and low-frequency regime 
[see Eq. 



ized by 

Abs = 4 - • (85) 

Michclson imperfections can be characterized by the dif- 
ference in the phase shifts and losses when light travels 
from the beamsplitter to the input mirrors of the two 
arms: 

ow.b = a M ± -Aa u , e MA B = e M ± ■ (86) 

In addition to A BS , Aa M and Ae M , which concern the 
beamsplitter, we consider the following contributions to 
mismatch between the arms, 

T. lA . B = Ti±^A T , (87) 

T eA , B = r e±^A £ , (88) 

<Aa,b = 0±^A , (89) 

that is, mismatch between input mirror power transmis- 
sivities, end mirror losses and cavity detuning, respec- 
tively. We replace these with the following more conve- 
nient quantities: 

Ae A T Ae L = A e AA _ A^ 
e Ti +T e ' e Ti + T e ' X <j> ' [ ' 

[See Table ITT1 for definitions of e, e L and A.] 

In the remainder of this section, we give the trans- 
fer functions from the carrier light (DC), laser amplitude 
fluctuations and laser phase fluctuations to the differ- 
ential output port, to first order in the mismatch (re- 
call that ideally, in the absence of imperfections, these 
common-mode inputs do not appear in the differential 
output port). We keep our formulae to the leading order 



in {QL/c, eL/c, XL/c}, and ignore the averaged losses, e L 
and e M (but not Ae L and Ae M ). We refer to this as the 
leading-order approximation. Furthermore, in order to 
keep the analytical results understandable, we work only 
in the regime of {57, 9 C } -C {0 D ,A, e}, which we shall 
refer to as the low-frequency regime. 

Definitions and assumed values for A BS , Aa M , Ae M , 
A T , A<£, and A e are given in Table ||| 



1. The Carrier 

The transfer function from the carrier to the differen- 
tial output can be written as 

x>^(::|). m 

where definitions of A( fc ), values of and y>£L are 
listed in TablellVl assuming the carrier at the beamsplit- 
ter is in the first (amplitude) quadrature. 

Contributions listed in Table IIVI can all be obtained 
from simple considerations. First, since each field that in- 
terferes at the beamsplitter is scaled by one transmission 
and one reflection coefficient factor, A BS does not con- 
tribute to the output carrier light at the differential port. 
Then, for all mismatches except the loss, one only has to 
notice that when the arm cavities are lossless, carrier light 
with amplitude D and phase ip = returns to the beam- 
splitter with amplitude reduced to (1 — e M ), and quadra- 
ture rotated by 2a + 2a M . As a consequence, the differen- 
tial output port gets (D/2)(— Ae M ) = (— Ae M /2)D in the 
tp = 2a quadrature (factor of 2 due to the beamsplitter), 
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and 



(D/2)A[2a + 2a M ] 
eX ( Ae AX 

T + ~ 



x 2 



Am 



D 



(92) 



in the orthogonal quadrature, ip = 2a+n/2. The effect of 
the loss mismatch can be understood when we decompose 
the (complex) reflectivity of the cavity into a sum of two 
components: 



*7Le 2ict > 



'Rj 



1 - ^/RjTee 2 ^ 



1 + iX/e e L 2 
1 — iX/e el— iX/e 
2ia £ l 2e ia 



(93) 



Here we see that the loss e L creates an output at the ip = 
a quadrature, so an imbalance in loss Ae L will contribute 



Ae T 



D 



(94) 



e Ve 2 + X 2 

in the tp = a quadrature in the differential output port. 

2. Amplitude (Intensity) and Phase (frequency) Noise 

Under our simplifications, the laser amplitude noise 
Z\ and phase noise z 2 couple to single (yet frequency 
dependent) quadratures in the differential output port, 
as parametrized by 



k 



N, 



(fc) 



(k) 



- sin pf k} 
cosp( k) 

-sin^f fc) 
costf k) 



Zl 



Z2 



(95) 



Measurement of the output quadrature = b\ cos £ 
&2 sin£ will include the laser noise: 



Nfczi sin(vfk) ~0 + Nj k) z 2 sin(v>f fc) - C) 



(96) 

In particular, the output quadrature C, = is not sen- 

sitive to the fc-th contribution of laser amplitude (phase) 
noise [note that we have switched the notation for ip from 
that of Eq. @3J]. 

As it also turns out, in the leading-order approxi- 
mation and the low-frequency regime, NEycpf^ = 

\_C(k) ' f( k)] ■ Considering the different ways tp appears 
in Eqs. (j91(l and l|95|) . this means the phase noise cou- 
pled to the differential output port remains orthogonal 
to the carrier. This can be argued for easily: since phase 
modulations on the carrier do not drive mirror motion, 
the propagation of phase noise is not affected by the op- 
tical spring. Amplitude modulations, on the other hand, 



do drive mirror motion and therefore should couple to 
the differential port in a dramatically different way. We 
tabulate the quantities and ip^ in Tab. from 

which we can see that the amplitude-noise coupling has 
features around the common-mode optical-spring reso- 
nant frequency, <d c . 



3. Evading Laser Noise by Artificial Asymmetry 

For realistically achievable symmetry between the two 
arms, laser noises turn out to be the dominant noise 
source to our squeezer. Here we discuss a novel way 
of mitigating laser noise coupling by introducing artifi- 
cial asymmetries. According to the approximate results 
(in the leading-order approximation and low-frequency 
regime) obtained in the previous section, both ampli- 
tude and phase noise emerge from single quadratures (as 
vector sums of contributions from different mechanisms). 
We can, therefore, eliminate the laser noise totally, up to 
this order, if we make both of them emerge from the same 
quadrature £ + n/2, and make sure that the orthogonal 
quadrature, £, has a sub- vacuum noise spectrum. At our 
disposal are two asymmetries that we can adjust manu- 
ally: Aa M and Ae M . 

At any given sideband frequency 17, for a generic set 
of other asymmetries, it is always possible to make both 
laser noise sources emerge at the C+^/2 quadrature (and, 
therefore, to vanish at the £ quadrature), by adjusting 
Aa u and Ae M , if the following non- degeneracy condition 
is satisfied: 



det 



r(fi,C) 



sin (¥>L M 



* 0. 



(97) 



[See Eq. Qfy.] 

According to Tables IIVI and laser phase noise 
emerges in a frequency-independent quadrature, but the 
amplitude noise does not. This means the elimination 
of laser noise must be frequency-dependent, and we can 
only choose one particular frequency for perfect laser 
noise evasion. However, if £1 ^> C is also satisfied, 
then the frequency-dependence goes away. We consider 
this special case, and choose a detection quadrature of 
C = 3a/2, i.e., the one with minimum quantum noise. 
From Tabs. HVl and Ivl we get 



' 2 



X 2 



?0. 



Since the carrier always emerges n/2 away from the phase 
noise, it emerges in exactly the same quadrature we 
propose to detect. In this way, the laser- noise-evading 
squeezer always produces squeezed light with amplitude 
squeezing. 
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Finally, we note that, due to possible higher-order cor- 
rections, laser noise evasion may not be as perfect as pre- 
dicted by our first-order approximation, even at a single 
frequency. The amount of residual laser noise, as well as 
the exact level of the deliberate asymmetries we intro- 
duce, must be given by a more accurate calculation. 



C. Comparison between analytical calculations and 
numerical simulations 

In Table [IJ we list the parameters used in modelling 
our interferometer. An important feature of the numeri- 
cal code is that it can handle imperfections in the optics 
quite naturally, while for analytical techniques the solu- 
tion becomes complicated rather dramatically when more 
ingredients are added. To fully test this feature, we con- 
structed a test case with realistic imperfections. The im- 
perfections included were those mentioned in Sec. IVI Bl 
Using the parameters listed in Table [I] we calculate the 
noise at the differential port due to quantum fluctuations 
entering from this port and from lossy mirrors, as well as 
laser amplitude and phase fluctuations entering from the 
symmetric port. 

In Fig. we show the calculated noise levels from nu- 
merical simulations in curves, while those from the ana- 
lytical treatment are shown as solid points. The agree- 
ment between the two sets of calculations is reassuring. 
Now we discuss these noise spectrum in more details. 
In the left panel of the Fig. [7| we plot noises due to 
vacuum fluctuations entering from the dark port (light 
curve and points), and due to vacuum fluctuations en- 
tering from mirror losses (dark curve and points). In 
both results, there is a rather dramatic resonant fea- 
ture around the differential-mode optical-spring resonant 
frequency, at D ~ 8 kHz, as can be expected from 
Sec. IVI Al The rather weak but still noticeable fea- 
ture around the common-mode optical-spring resonant 
frequency O c ~ 360 Hz is solely due to optical parameter 
mismatch. In the right panel, we show laser amplitude 
(light curve and dots) and phase (dark curve and points) 
noises; we have introduced artificial asymmetries a M and 
e M , with values obtained empirically using the numeri- 
cal simulation code, such that both laser noise sources 
are evaded to a roughly maximal extent at 1 kHz. For 
this reason, contributions to the results shown here are 
largely higher order, and we cannot hope to explain them 
using results obtained in Sec. IVIBI Here we do observe 
dramatic features around both the differential-mode and 
the common-mode optical-spring resonances. 

Results in Fig.[7|are also of great significance for a prac- 
tical reason: they show that the vacuum modes exiting 
the interferometer are squeezed by a large factor even in 
the presence of realistic estimates for optical losses (left 
panel) and laser amplitude and phase noise (right panel) . 



VII. SUMMARY AND CONCLUSIONS 

The main purpose of this work was to develop a mathe- 
matical framework for the simulation of quantum fields in 
a complex interferometer that includes radiation pressure 
effects. We work in the linear regime around the oper- 
ation point of this interferometer; in this regime, after 
adopting the Heisenberg picture of quantum mechanics, 
the quantum equations of motion (Heisenberg operators) 
of observables are identical to classical ones. 

During the development of this framework, we aug- 
mented previous treatments of mirrors (and beamsplit- 
ters) by allowing the carrier phases at the four (eight) 
ports to be different. This extension gives rise to the op- 
tical spring effect even without detuned optical cavities. 

Based on this mathematical framework, we developed 
a simulation code that can allow arbitrary optical topolo- 
gies, and applied it to a specific example of the interfer- 
ometer shown in Fig. [3] This interferometer was shown 
to be capable of squeezing the vacuum modes that en- 
ter - and subsequently exit - the differential port of the 
beamsplitter. We introduced optical spring effects by 
detuning the arm cavities as a means of mitigating the 
detrimental effects of thermal noise. We study not only 
the quantum noise, but also laser noise couplings from 
the symmetric (input or bright) port to the output (an- 
tisymmetric or dark) port. Good agreement was found 
between numerical results given by this code and analyt- 
ical ones derived independently. This agreement makes 
us confident that the simulation is working correctly for 
this rather complex interferometer. 

During our study of the laser noise couplings, we found 
a novel method of evading the laser noise by introducing 
artificial but controlled asymmetries. This is crucial for 
the practical implementation of this interferometer, and 
is likely to find applications in many other experiments. 

Our simulation code is now being used in the detailed 
optical design of the Advanced LIGO interferometer. We 
also envisage the following extensions to the code in the 
near future: 

• Allowing multiple carrier or rf sidebands, which 
may be relevant to the modeling of squeezing ex- 
periments that use nonlinear optical media, e.g., 
crystals, as well as the modeling of error signals for 
control systems. 

• Incorporating the modeling of servo loops. Here we 
may rely on the input from quantum control theory 
as to whether and how realistically a Heisenberg 
treatment can describe a electro-optical feedback 
system. 

• Allowing nonlinear media or other elements with 
"custom" dispersion relations. 
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FIG. 7: Spectra of noise power at the output port of the ponderomotive interferometer, normalized to the vacuum noise level. 
The noise power is dimensionless, as compared to vacuum; a pure vacuum (or shot noise) corresponds to unity. The lines are 
results of the simulation code, while the data points are values calculated from the corresponding analytic calculations. In 
the left panel the vacuum noise level of the light exiting the antisymmetric port is shown. The solid (blue) curve shows the 
vacuum noise due to the (unsqueezed) vacuum fluctuations that enter via the antisymmetric port of the interferometer; the 
dashed (green) plot represents the noise due to the vacuum fluctuations that enter via other optical losses in the system. At 
all frequencies where the vacuum noise power is below unity, the vacuum modes exiting the interferometer are squeezed due to 
radiation-pressure effects. For the squeezing to be useful, all noise couplings must yield a lower noise power than the squeezed 
vacuum. In the right panel we show the coupling of laser frequency (solid, blue) and laser amplitude (dashed, green) noise fields 
to the output port, as calculated by the simulation code. Noise levels of 10 -4 Hz/Hz 1//2 for frequency noise and 10 _8 Hz -1 ^ 2 
for amplitude noise are assumed at the input to the interferometer; all other parameters are listed in Table [I] 
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APPENDIX A: TWO-PHOTON QUANTUM 
OPTICAL FORMALISM 



where Ei(t) and E 2 (t) are called quadrature fields, which 
vary at timescales much longer than that of the optical 
oscillation, 1/cj. The quadrature formalism replaces E(t) 
by 



E 



Ei{t) 
E 2 (t) 



(A2) 



The DC components of Ei >2 (t) can be regarded as 
monochromatic carrier light. In particular, carrier light 
with amplitude De lip is represented as 



De l(p (De lv ) e~ iut O D 



COS if 

sin tp 



(A3) 



AC components of Ei j2 (t), which we denote by Ai i2 (t), 
are called sideband fields, which are usually more conve- 
nient to study once transformed into the frequency do- 



We use the two-photon formalism developed by Caves 
and Schumaker [Hj to describe GW interferometers 
with significant radiation-pressure effects. In this for- 
malism, any quasi-monochromatic optical field A with 
frequency near the carrier frequency lu is written as 



E(t) = Ei(t) cos(wt) + E 2 (t) sm(wt) 



cos tot sinwi 



(Al) 



ii, 2 (fi) 



A h2 (t)e mt dt 



(A4) 



In quantum two-photon optics, it is convenient to use 
a particular normalization for sideband fields: 



Ac 



[ +C °^[a h2 (n) e -^ + H.c.] 
Jo Zn 



(A5) 
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In this way, we have a convenient set of commutation Here we have denoted a\ i2 = ai^iS^), a'\ 2 = &i.2(^')- 
relations (for Q < lj) |l[Tol|: 



[oi, a[) = [a 2 ,a' 2 ] = [01, a'{] = [a 2 ,a^] = , (A6a) 
[01, 4] = -[03, ai f ] = 27ri(5(Q - fi') . (A6b) 
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